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We consider deuterium compressed to higher than atomic, but lower than nuclear densities. At 
such densities deuterium is a superconducting quantum liquid. Generically, two superconducting 
phases compete, a "ferromagnetic" and a "nematic" one. We provide a power counting argument 
suggesting that the dominant interactions in the deuteron liquid are perturbative (but screened) 
Coulomb interactions. At very high densities the ground state is determined by very small nuclear 
interaction effects that probably favor the ferromagnetic phase. At lower densities the symmetry of 
the theory is effectively enhanced to 5*17(3), and the quantum liquid enters a novel phase, neither 
ferromagnetic nor nematic. Our results can serve as a starting point for investigations of the phase 
dynamics of deuteron liquids, as well as exploration of the stability and dynamics of the rich variety 
of topological objects that may occur in phases of the deuteron quantum liquid, which range from 
Alice strings to spin skyrmions to Z2 vortices. 



I. INTRODUCTION 

It was recently pointed out PQ that deuterium, when 
compressed between atomic and nuclear densities, be- 
comes a quantum liquid over a large range of temper- 
atures. This comes about through the confluence of a 
number of factors. First, when the interparticle dis- 
tance I is smaller than twice the Bohr radius oq, deu- 
terium atoms overlap and are ionized. One then has 
two separate fluids, one composed of electrons and an- 
other composed of deuterons. Due to the Coulomb re- 
pulsion between deuterons, which is screened by the 
electrons only over distances much larger than /, the 
deuterons at zero temperature will crystallize into a lat- 
tice. The temperature at which the crystal melts scales as 
o-oTcrys ^ 1/180 XQ;(ao//), where a is the electromagnetic 
coupling constant. On the other hand, the deuterons 
Bose-condense at temperatures below aoTcond ~ 47r^/3 x 
{Mao)~^{ao/l)'^, where M is the deuteron mass [5]. (We 
work with units where ks = fi = c = 1.) Therefore, at 
high enough densities {I < ag), there is a range of tem- 
peratures, Tcrys < T < Tcondi wherc a quantum liquid 
of deuterons should form. Of course, it must also be the 
case that the relevant densities are still far from nuclear 
ones, and Ref. [T] pointed out that this is indeed the 
case for deuterium. Since deuterons are charged bosons, 
the quantum liquid will be a superconducting superfluid. 
This is the regime we address in this paper. 

There are two motivations for studying this system. 
First, this kind of matter is expected to exist in a layer 
inside brown dwarfs that are light enough not to ignite 
deuteron fusion [1. Second, it may be created in terres- 
trial laboratories through shock compression [31 13], ex- 



periments using inertial confinement or other techniques 
[S1[B]. To get an estimate of the pressures required, sup- 
pose that the pressure is dominated by the electron de- 
generacy pressure. This should be reasonable once the 
deuterium is ionized, which one would expect to take 
place once / ^ 2ao. The pressure P is then given by 



P w 1.8 X 10' 



,3 /2ao\5 



GPa, where we define the interpar- 
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tide distance as I = m}^^, and n is the particle density. 
Since Tcond exceeds Tcryst once 2ao/Z ~ 1, deuterium 
should become a quantum liquid of deuterons once the 
pressure reaches Pg ^ 1.8 x lO'^ GPa. Pq sensitively de- 
pends on the estimates for the ionization density and 
Tcond, Tcryst, SO that Pq should be viewed as only a rough 
estimate. For instance, recent quantum molecular dy- 
namics simulations suggest that deuterium may already 
be ionized at densities of ^ 400 GPa [7] , an order of mag- 
nitude less than naively setting I = 2ao above would sug- 
gest. 

Current diamond anvil cells reach a pressure of a few 
hundred GPa, as do inertial confinement experiments [51- 
IT2] . Both are a factor of 10— 100 shy from our estimate of 
the pressures necessary for the appearance of the phases 
which we will be discussing. However due to the im- 
portance of the subject for energy production and other 
applications, one can hope that this gap will be closed in 
the future. Our goal is to explore some of the unusual 
properties of this kind of matter as a guide to these ex- 
periments. 

Since we are interested in phenomena occurring at 
length scales much larger than the size of the deuterons, 
the system can be described by a Lagrangian with point- 
like deuterons and electrons interacting via electromag- 
netism. In the regime of interest, the momenta of 
the deuterons are very small compared to their mass 
M ~ 2 GeV, so the deuterons are non-relativistic. Since 
the deuterons are charged spin-1 bosons, they can be 
described by a complex 3- vector field d. The non- 



relativistic Lagrangian governing this system is 
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C + tiN ^d^ {iDq + /i + — ) d + V(J^+ m + /ie7°)^ 
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B • (d^ X d) 



(1.1) 



where 13^ = d^ + iev4^, and ip and A^ are, respec- 
tively, the electron and photon fields, B is the magnetic 
field, g ~ 0.857 the magnetic moment of the deuteron in 
Bohr magnetons, m is the electron mass, and /x, /ie are 
the chemical potentials of the deuterons and electrons; 
d^ X d = e^^''d*dk in a slight abuse of notation. The 
Lagrangian has an 0(3) global rotation symmetry, and a 
U{1) electromagnetic gauge svmmetrv|47|. Not shown 
explicitly in the Lagrangian are the deuteron interac- 
tion terms, and terms with more derivatives, the form of 
which is constrained by the requirement of Galilean in- 
variance. The inclusion of interactions among deuterons 
is complicated, since two deuterons at rest lie above the 
threshold for "^i/e -I- n,^H + p production, so that all 
of these different nuclei need to be included in an ef- 
fective field theory. This can be done following the ap- 
proach to such effective theories described, for instance, 
in Refs. [T31 [H] ■ These contributions are subleading but 
will, in fact, play a decisive role below. We postpone a 
further discussion of the deuteron-deuteron interactions 
until Sec. Uni 



II. PHASES AND THE EFFECTIVE 
POTENTIAL 

At low temperatures (T < Tcond ~ 10^ K for I ss lOao), 
the deuterons condense, leading to a ground state expec- 
tation value for the d field. We assume that the conden- 
sate is spatially homogeneous, and leave the analysis of 
non-homogenous phases for a later publication. The ex- 
pectation value of d can be split into real and imaginary 
parts as (d) = a-|- ib, so that a and b are real 3- vectors. 
Depending on the relative orientation of a and b, one 
finds two different phases, sometimes named 'ferromag- 
netic' and 'nematic' phases in the condensed matter/cold 
atoms literature. If the condensate is such that a and b 
are parallel, we have the nematic phase 

(d)=ae*". (2.1) 

In the nematic phase, a rotation of a by tt can be undone 
by sending a — ?> a + tt, so the order parameter is a 'di- 
rector' rather than a vector. If instead a and b point in 
different directions, we have the ferromagnetic phase 



(d) = a + ib, a X b 7^ 0, 



(2.2) 



where the spin S — ^d^^ xd = axbof the condensate 
is non- vanishing. 

Our first task will be to decide which of these two 
phases is actually realized in the deuteron liquid. To do 



FIG. 1: Example of two graphs whose derivatives evaluated 
at dM = n cancel in the effective potential. The dotted lines 
are electrons, the solid lines deuterons and the wavy line an 
electron-dressed Coulomb photon line. 



this, one must evaluate the effective potential y(d^,d) 
and minimize it. (y(d',d) should be seen as a function 
of the expectation values (d) and (d'''); to simplify no- 
tation we will drop the brackets except where this could 
cause confusion, but it is important to keep in mind that 
the arguments of F(d^^,d) are classical, not quantum, 
fields.) 

The effective potential F(d''', d) is given by the sum of 
all one-deuteron irreducible diagrams with d, d^ external 
legs and zero external momentum|15| . Due to rotation 
and U{1)em gauge symmetry it can only depend on the 
combinations (d^' • d) — a^ + b'^ and (d^' • d'f)(d • d) ::= 
(a^ — 5^)^ -|- 4(a • b)^. If we retain only the terms up to 
quartic order in the effective potential we have: 



y(d^d) 



-/idM + ci(d^d)2 -I- cad^^d^ 



(2.3) 



For C2 > 0, V^(d^^, d) is minimized by a = &, a-b = 0, that 
is, we have the ferromagnetic phase. In the opposite case, 
C2 < 0, l^(d''',d) is minimized by a parallel to b, which 
yields the nematic phase. Taking into account terms of 
higher order in d ' , d other intermediate angles between 
a and b may be favored. 

We will now argue that, at leading order in the expan- 
sion parameter l/ao (ao = l/arrie is the Bohr radius), the 
effective potential is given by an infinite series of one- 
loop diagrams containing only screened Coulomb pho- 
tons. Since the dominant Coulomb interaction is spin- 
independent, this amounts to the claim that ci 3> C2 



in Eq. (2.3). To show this we have to argue that all 



other diagrams are suppressed by powers of l/ao or by 
the deuteron mass AI. What complicates this argument 
is that the estimate of any diagram we neglect depends 
on the value of the chemical potential fi. However, re- 
call that fi serves to enforce charge neutrality, since it is 
is the chemical potential necessary for the existence of 
the same density of deuterons and electrons. This means 
that one can only compute /i after ^(d^^, d) is known. We 
break this impasse by positing a certain estimate for /i, 
and then making sure it is self-consistent. Specifically, we 
suppose that fi ^ a/ 1, and compute y(d',d) assuming 
this estimate of /i. We then use the calculated y(d^^,d) 
to compute /z by demanding that charge neutrality be 



enforced, which is the condition that 
9F(dt,d) 



9(dtd) 



= 0, 



dtd=ri 



(2.4) 



where n = l/l^ is the electron density, and check for 
self-consistency. 

Our first step is to integrate out the Coulomb photons 
Aq (ignoring the magnetic photons for the moment). The 
result is 






(2.5) 



where Vc is the Coulomb potential screened by the pres- 
ence of the electron Fermi sea; the Fourier transform of 
Vc{q) at small momenta [q ^ l/l) is 



VM 



(2.6) 



with the Debye screening mass m^ ^ am^/l = I/qqI [IB] . 
We do not include magnetic interactions at this stage 
since they are suppressed by powers of 1/M, and thus one 
would expect their contribution to the effective potential 
to be suppressed. Of course, this needs to be checked 
self-consistently, and the effects of magnetic interactions 
are discussed in Sec. IIIII 

Eq. (2.4 1 implies a cancellation between certain dia- 
grams. Consider the two diagrams in Fig. [T] The first 
one contributes to the dM term in y(d^^,d), while the 
second contributes to the (dM)^ term. Their derivatives 
evaluated at d^d = n, however, cancel against each other: 



d 



(9(dtd) 



;dtd 



a 
2m\ 



(dM)^ 



(2.7) 



The same cancellation occurs for any graph containing a 
"dangling" deuteron line connected to the rest of the dia- 
gram by only one photon line, and thus we can disregard 
them from now on. 

This cancellation has a simple physical interpretation. 
Each deuteron interacts with the average charge of the 
background of other particles. The charge of this back- 
ground, however, vanishes as the electron charge cancels 
that of the deuterons. Any shift in the energy will come 
from charge density fluctuations, and these are described 
by two (or more) photon exchanges. An example of a 
contribution of this kind is shown in the first graph of 
Fig. [2j More complicated contributions involving more 
deuteron and photon lines, however, contribute equally 
and need to be resummed. 

Fortunately, the resummation involved is fairly stan- 
dard, and amounts to the computation of the one-loop 
effective potential [15 . Physically, it can be seen as the 
sum of zero-point energies of the deuteron quasiparticles. 




+ / d3r'(dTd(r)-V^7>(r))Vc(r-r')(dTd(/)-?/'7V(?-')) 



FIG. 2: Class of graphs contributing to the leading order 
effective potential displayed in two different ways. 



including the Coulomb interaction with the condensate. 
Alternatively, it can be seen as the sum of the zero-point 
energies of the photon quasi-particles dressed by inter- 
actions with the electrons and the deuteron condensate. 
The one-loop effective potential Vi can be obtained by 
writing the deuteron field as d = (d) + Xi keeping only 
the terms quadratic in the fluctuations x^ , x and per- 
forming the resulting functional Gaussian integral. The 
result has the well known "trace log" form that gives 



5yi(dt,d) 



9(dtd) 



-M 



2 /■ (Pq K(g)(nye(g) + 2e,) 

3 7 (2^)3 E, 

(2.8) 
where e, = q^/2M - fi, E^ = e^ + {A/3)nVciq)eq. The 
electron density n = \/l^ appears here because, up to 
terms suppressed by higher powers of a, we can make 
the substitution ■ip^'^il; — >• n in Eq. (2.5 1. Note that the 



expansion of Vi in powers of Vc{q) gives the sequence of 
diagrams portrayed in Fig. [21 



The integral appearing in Eq. (2.8) would be infrared 
divergent if not for the electron screening (the mj. term) 
and condensate screening, which appears as the term pro- 
portional to nVc{q) in the denominator. That means that 
the integral will be dominated by small values of q, on 
the order of q ^ m^. When I <C ao and /i ^ a/l we 



have that £„ 



M 



a/l ~ {l/ao)na/mj ^ nVc, and so 



eq can be neglected in the numerator. Similarly, in the 
denominator, we use the hierarchy of scales 

q'^/M ~ m^jM < a// - ^ < na/m^ - nVc- (2.9) 

This means that the condensate-screening nV{q) term in 
Eq in the denominator gives the dominant contribution 
in the regime of interest. We see then that the infrared 
divergences in the potential are partially cutoff by the 
Debye screening mass m^, and partially by the effect of 
the condensate. 

Using the results above, we can now estimate the in- 
tegral in Eq. (2.8) and find the value of the chemical 



potential necessary to keep the deuteron density equal to 



the electron density. This condition turns out to imply 
that 



M' 



nor Ivfi^g 1 
y//Z V na nis 



M' 
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7' 



(2.10) 



so that our estimate of /i is self-consistent, as previously 
advertised. 

Notice that each individual diagram that was summed 
above is parametrically larger than their sum. In fact, 
the one-loop diagram with n photon lines is of order 
(a/^'*)(ao/Z)"~'^/^ so they are larger the more photon 
lines they have. Their sum, however, is smaller, of order 
a//*. This is because the infrared divergences are actu- 
ally cut off by the scale set by the deuteron condensate, 
an effect not included at any finite order in perturba- 
tion theory, instead of the scale /i appearing in individual 
Feynman diagrams. 

What remains to be done now is to argue that other 
graphs would, in the regime considered, give contribu- 
tions smaller than the ones that we kept. As mentioned 
before, these estimates hinge on the value of ^ ~ a/l and 
finding this value was indeed the main motivation for the 
effective potential argument above. The suppression of 
the remaining diagrams is best argued through some ex- 
amples. Consider, for instance, adding one photon line to 
the one-loop diagrams of Fig. [2j The effect is to substi- 
tute one Vc{q) by a two-photon ladder with zero incoming 
(relative) momentum and q ~ ms outgoing momentum. 
This ladder diagram can be estimated as 



d-^k 



1 



g2 -I- to2 (^q _ ]^Y + 



ml efc 



^inis 




(2.11) 

and, consequently, gives a smaller contribution compared 
to the leading diagrams by a factor of ^JiJoq. Graphs 
including electrons are also suppressed. For instance, let 
us look at at graph in Fig. [3] Its contribution to the 
effective potential can be estimated as 



(2.12) 



where we used the fact that the electron loop, at small 
momenta q ~ rus is of order ~ mil |16j . This contri- 
bution is again suppressed by yjl/ao compared to the 
leading order. The picture emerging then is very simi- 
lar to the high density limit of the jellium model, which 
is a charged Bose gas with a fixed, non-dynamical back- 
ground of negative charges [T71 [TH] , modified by the elec- 
tron screening of the Coulomb force. 



III. SYMMETRIES OF THE DEUTERON 
LIQUID 

The most important consequence of identifying the 
one-loop Coulomb diagrams as giving the leading con- 





FIG. 3: An example of a sub-leading graph. The dotted line 
represents an electron. 



FIG. 4: Leading electromagnetic SU{2>) breaking effect. The 
square vertex represents a magnetic dipole interaction. 



tribution to the effective potential is that the effective 
potential is a function of d^d only, not of d^ d^. This 
can be seen from the fact that all the incoming deuteron 
lines are contracted with outgoing deuteron lines on the 
graphs in Fig. [2] Physically, this follows from the spin- 
independence of the Coulomb forces. This implies that 
the symmetry group of our theory, which is 0(3) x 
U{1)em, is effectively enlarged to SU{Z) x U{1)em HE]. 
We now examine the size of the SU{3) breaking effects. 
The leading electromagnetic SU{3) breaking diagram is 
shown in Fig. [4] This diagram can be estimated as 
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id*d,d' 



didjd'' - 2d*djd^ • d) 



(3.1) 



The loop integral above, which would be IR-divergent 
if not for static and dynamic screening effects, con- 
tains three scales: /i, ttIj, and the scale associated with 
the Landau damping of the magnetic photon, q^ ^ 
{m^l^y/'^. For q ^ qm, the propagator of the magnetic 
photon becomes qtqj/q'^ — >■ Sij. For smaller g, qm sets 
the scale for the momenta in the magnetic photon prop- 
agator. Using the fact that fj, ~ a/l < ms ~ (aoZ)^"^/^, 
the estimate of the integral is 
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break 
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(3.2) 



This contribution would favor the ferromagnetic phase 
over the nematic phase. Notice that the integration over 



q sets q^q^ /q^ — >■ (J'-'/S, implying that the only part of 
the dipole potential surviving is, in position space, a 
5{r) function at the origin. This leads to the surpris- 
ing conclusion that contact interactions due to nuclear 
forces may compete with dipole-dipole interactions be- 
tween the deuterons. A very similar situation involving 
dipole-dipole interactions occurs in the physics of W bo- 
son condensation in the early Universe [H] . There the W 
boson contact interactions are fixed by gauge symmetry 
and are known to favor the ferromagnetic phase. 

In our case the the strong nuclear interaction between 
two deuterons is more uncertain. In the effective theory 
we are using, valid at distances larger than the deuteron 
size (w 5 fm), the s-wave deuteron-deuteron interaction 
is described, at lowest order in the momentum expansion 
by 



-\7-nuc 
^ break 
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(3(dt.d)2 



- dt'd' 



(3.3) 

The term proportional to Aq {A2) contributes to the spin 
channel S — Q [S = 2) deuteron-deuteron channel. In 
fact, in our effective theory, the s-wave scattering lengths 
in the spin S — 0,2 channels would be given by the terms 



in Eq. (3.3) dressed by Coulomb photons with momenta 
within the range of validity of the theory [Q < 1/5 fm~^), 
which are explicitly included in Eq. ( 1.1 ). Thus, the coef- 



ficients Aq,A2 above are the deuteron-deuteron scatter- 
ing lengths in the absence of soft Coulomb interactions, 
but including the effect of hard photons {Q > 1/5 fm~^). 

What complicates the description of deuteron- 
deuteron scattering is the existence of channels below the 
deuteron-deuteron threshold, namely, n+^He and p+'^ H . 
One way of dealing with these additional channels is to 
include n,p,^He and ^H fields explicitly in the effective 
theory, together with their respective couplings, and per- 
form a coupled-channels calculation in the medium. Most 
likely this is unnecessary. The presence of these open 
channels will make the deuteron-deuteron phase shifts 
complex, but this can be rigorously taken into account 
by taking the parameters Aq,A2 to be complex. This 
technique is used, for instance, in non-relativistic QCD 
where a quark-antiquark system can annihilate into en- 
ergetic gluons [20^. 

In any case, one expects Aq and A2 to have very small 
imaginary parts. To see this, note that the deuteron- 
deuteron initial state can have S* = or 5' = 2. First, 
consider the S — 2 deuteron-deuteron initial state. The 
intermediate state can have only S* = or S = 1, so 
it will have some non-vanishing angular momentum. At 
small energies, mixing with higher partial waves is sup- 
pressed, so the scattering length in the S = 2 channel 
should have only a small imaginary part. Next, consider 
the S = Q initial state. This initial state can mix with the 
intermediate S = L — Q state. However, the wave func- 
tion configuration of the initial and intermediate states 
is very different, so this mixing should be very small. 

In fact, a model calculation of the S" = deuteron- 



deuteron scattering length gives Aq = 4.91 ± 0.02 -|- 
i(-0.0115 ± 0.0001) fm [H]. In our effective theory, the 
values of Aq and A2 are the scattering lengths in the spin 
S' = 0, 2 channels in the absence of Coulomb interactions. 
More precisely, Aq and A2 subsume electromagnetic in- 
teractions at distances smaller than « 5 fm (or photons 
with momenta ^1/5 fm); the contribution from softer 
photons is included explicitly in Eq. (1.1). 



If the deuteron forces were spin-independent, Aq would 
equal A2, and we would have an effective potential with 
the enhanced S'C/(3) symmetry. Current few-body nu- 
clear technology is available for a realistic calculation of 
these parameters, but we are aware of only two published 
works. One model calculation, using a simple Malfliet- 
Tjon potential ^22j gives Aq — 10.2 fm and A2 — 7.5 
fm [23j. These values for the scattering lengths lead to a 
ferromagnetic phase. Ref. [H] computes only A^ giving 
the value quoted above {Aq = 4.91 ± 0.02 -I- i(-0.0115 ± 
0.0001) fm). 

A. Deltuva was kind enough to use the methodology 
described in [53] , |55] to compute the required scattering 
lengths without electromagnetic forces at our request |26j . 
The results he finds are Aq = 5.35,5.13,4.87 fm and 
A2 = 3.16,3.16,3.18 fm, respectively for the AV18, CD- 
Bonn and INOY04 potentials. The imaginary parts are 
all, as expected, of the order of 1%. None of these cal- 
culations include the effect of three-body forces, known 
to contribute a small (about 5%) to the binding energy 
of small nuclei or the hard photon exchange. There is a 
high degree of universality in low-energy few-nucleon re- 
actions - understood from the point of view of effective 
field theories — so the discrepancy of this result with the 
one in |53] can hardly be blamed on the difference be- 
tween models used. In any case, an estimate of the size 
of these scattering lengths can be made by just assum- 
ing a scattering length comparable to the "size" of the 
deuteron itself, namely, a few Fermi. This means that 
^break ~ A/M with A w 5 fm is actually larger than 
Kl^f^, since 



break'> 



M2 ^ M 



(3.4) 



given that 1/M -- 10^^ fm. 



It is natural to ask whether the nuclear interactions 
can favor angles between a and b that are not or 7r/2. 
As noted after Eq. (2.3), such an effect would have to 



come from terms with more than four powers of the d 
fields in the effective potential. Of course, the one-loop 
effective potential contains terms with all powers of d^ -d, 
and it is spin-independent. The leading spin-dependent 
contribution to the effective potential can be obtained by 
replacing one of the Coulomb photons in the diagrams 
leading to the one-loop effective potential with a nuclear 
contact interaction between the deuterons, as shown in 
Fig. [5] However, it is not hard to see that this merely 




FIG. 5: Leading nuclear 5'[/(3)-breaking contribution to the 
effective potential. The contact term is a spin-dependent nu- 
clear interaction. 



gives a contribution to the potential of the form 



'"^V(dt.d)dt^d^ . 4"^^ 



3M 



3Af 



(d^-d) (3(dt 



dt'd^ 



(3.5) 

for some functions / and g that depend only on dM. 
Thus at this order the effective potential would pick out 
a phase with either a x b = or a • b = 0. To get 
other angles between a and b, we must have at least 
two nuclear contact interactions in the effective potential. 
This sort of contribution is highly suppressed by powers 
of \/M. So to leading order, the nuclear interactions 
drive the dcuteron liquid into either the nematic phase, 
or the ferromagnetic phase with a • b = 0. 

Finally, one may wonder how it is that nuclear ef- 
fects can come to affect the deuteron liquid, since the 
strong interactions are short-ranged, and naively the 
deuterons are prevented from getting close to each other 
by Coulomb repulsion, which is only screened on dis- 
tances larger than oq. However, in the regime of inter- 
est the dynamics of the deuterons cannot be thought of 
in terms of collisions between classical particles. The 
deuterons are in a quantum condensate, so that their 
wave-functions overlap, and they are highly off-shell, 
making classical intuition about scattering difficult to ap- 
ply. One might have thought that electromagnetic dipole 
interactions would be long-ranged compared to the nu- 
clear interactions, but as it turns out they yield only con- 
tact interactions to leading order, due to angle-averaging 
and the spatial homogeneity of the condensate. Since 
the strength of the electromagnetic SU(3)-breaking con- 
tact interactions turns out to be much smaller than that 
of the nuclear ones, it turns out that nuclear interac- 
tions can play an important role in the dynamics of the 
deuteron liquid. 



IV. FINITE TEMPERATURE EFFECTS 

Since the size of the SU{i) breaking effects is sup- 
pressed by powers of 1/M, one may wonder whether 
they will actually have any physical consequence after 



the small, but finite, temperature effects are included. 
One way of estimating this is to compare the energy den- 
sity required to move the a and b fields from parallel to 
orthogonal directions 



c 47rA 1 2 n 



(4.1) 



to the thermal energy density available in the system 

(4.2) 
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When £t is greater than £g, one would expect that 
SU(3)-breaking effects get washed out. We then find that 
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(4.3) 



At the typical temperatures (of order IQ^K => aoT ^ 1) 
that we are interested in, where T is larger than the crys- 
tallization temperature but smaller than the condensa- 
tion temperature, the criterion above splits the interest- 
ing range of densities into two regions. 

At larger densities {I < 10~^ao) the S'C/(3)-breaking 
effects are important and determine the phase of the sys- 
tem. With Aq > A2 the ferromagnetic phase is favored. 
Once the global SU{3) symmetry is broken, the sym- 
metry of the system prior to deuteron condensation is 
0(3) X U{1)em- In the presence of a ferromagnetic con- 
densate (d) = a. + ih, which has a • b = and a — b, 
the symmetry breaks to to a global Uz-q(1) subgroup 
combining a rotation around the a x b 'z' axis and an 
opposite phase rotation. The vacuum manifold is then 
(0(3) X U{1)em)/U,^q{1) ~ 0(3). This pattern of sym- 
metry breaking also occurs in cold atom optical traps of 
spin-1 atoms [57]. 

If it turns out that A2 > Aq, the nematic phase would 
be favored instead. In this phase the 0(3) x U{1)em 
symmetry is broken to 0(2) k Z2. If the condensate 
points along the 2 direction, the 0(2) consists of rota- 
tions around the z-axis, as well 2D-parity transformations 
(x, y) — > (— x, y), while the Z2 is realized as a rotation of 
the condensate by tt around the x or y axes, followed mul- 
tiplication by e" € U{1)em- Thus the vacuum manifold 
is (0(3) X C/(1)£m)/(0(2) k Z2). 

At lower densities, which are more accessible to ex- 
periments and observation, £g/£T -C 1 and the SU{3)- 
breaking effects are relatively unimportant. In that case 
the symmetry of the theory is effectively enhanced from 
50(3) xUil)EM to SU{3) X U{1)em- Using SU{3) x C/(l) 
transformations it is possible to go from the ferromag- 
netic phase with a-b = to the nematic ax b = phase. 
In this regime, the deuteron liquid is in a novel phase, 
which we will refer to as the SU(3) phase; it is neither a 
ferromagnetic nor nematic phase. The vacuum manifold 
can be found by noting that by a choice of coordinates in 
the three-complex dimension space we can make (d) to be 
real and have only its ^-component non-vanishing. The 



symmetry breaking pattern is then SU{3) x U{1)em — > 
C/(2), where C/(2) is the subgroup of C/(3) that leaves 
(d) invariant. The vacuum manifold is then a five- 
dimensional sphere {SUi3) x U{1)em)/U{2) = S^. 



V. TOPOLOGICAL DEFECTS 

The vacuum manifolds of theories with charged con- 
densates generally allow for the existence of topological 
defects. For the condensation of spin-0 nuclei, this has 
been examined in Ref. [55], following some foundational 
work in Refs. [H [5UH34) . The vacuum manifolds of the 
deuteron liquid, in which spin-1 nuclei are condensed, 
support a rich variety of finite-energy topological defects. 
We leave a detailed study of their stability, dynamics, and 
physical implications to future work, and confine our- 
selves here to simply sketching how they appear. 

Let us first discuss the situation for the more acces- 
sible densities at which Eb/Et <C 1, so that SU{?>)- 
breaking effects are negligible. The ground state mani- 
fold in this case is simply {SU{i) x U{1)em)/U{2) = S'^ 
which has trivial homotopy groups Hk until fc = 5. Since 
■"■5 ('5'^) = Z, the effective theory of the deuteron liquid 
should admit a Wess-Zumino-Witten term with a quan- 
tized coefficient [3S] . It would be interesting to work the 
physical implications of such a term in this context. 

It is sometimes assumed that finite-energy topologi- 
cal defects are classified simply by the homotopy groups 
of the coset space G/H, where G is the full symme- 
try group of the theory and H is the subgroup of G 
that leaves the condensate unchanged. For instance, one 
might think that if Tri{G/H) = the theory would not 
support vortex strings. However, the situation is actu- 
ally more subtle if only a subgroup Gl of G is gauged 
and breaks to a gauged subgroup Hl of H. When this is 
the case the homotopy groups of the coset space Gl/H^ 
of gauge symmetries become important. When these ho- 
motopy groups are non-trivial, various 'semilocal' defects 
are possible [36MT] . 

In our case, Hj^ is trivial and Gl = U{1)em, so 
the coset space is simply U{1), and the theory supports 
semilocal vortex strings because 7ri(t/(l)) = Z. To see 
why a non-trivial Gl/Hl can result in vortex strings, 
note that the action of the U{1)em on (d) defines a cir- 
cle in the vacuum manifold S^ . Motion along this circle 
costs no energy, and the value of the deuteron field at 
spatial infinity can wrap the circle in S^ defined by the 
gauge orbit, giving a topologically non-trivial configura- 
tion. The semilocal vortex string cannot unwind even 
though 7ri(S'^) = 0, because to do so it would have to 
leave the gauge orbit at infinity, which would cost an in- 
finite amount of gradient energy. This discussion makes 
clear that the non-trivial embedding of the gauge vacuum 
manifold Gl/ He = U{1) ~ S^ in the full vacuum mani- 
fold S^ gives rise to a fibration of S^ over CP^ with fiber 
S^. Since 7r2(Cp2) = Z, the S'[/(3) phase of the deuteron 
liquid can contain global monopoles; these monopoles live 



on the ends of semilocal vortex strings [351 135] • Semilocal 
defects are not always stable, and we leave an investiga- 
tion of their stability and physical effects using the power 
counting developed in this paper to later work. 

Let us now consider the two possible phases of the 
deuteron liquid at higher densities, I < lO^^ao), when 
the 5'[/(3)-breaking effects are important and determine 
the phase of the system. First, consider the ferromag- 
netic phase. The symmetry 0(3) x U{1)em is broken 
down to a global U(l) symmetry; as noted previously the 
global U(l) is generated by a linear combination of gener- 
ators of the global 0(3) symmetry and the local U{1)em 
symmetry. Before discussing semilocal defects, consider 
the usual kind of topological defects, which can be clas- 
sified by the homotopy groups of the full vacuum man- 
ifold. The vacuum manifold in the ferromagnetic phase 
is just 0(3). Since 7ri(0(3)) — Z2, the ferromagnetic 
phase supports topologically stable vortices, but ones 
that are rather unusual: rather than one vortex for each 
integer value of the circulation, there is only one kind 
of topologically non-trivial vortex. Furthermore, since 
713(0(3)) = Z, the ferromagnetic phase supports topo- 
logically stable 'spin Skyrmions'. (Adapting the usual 
hedgehog ansatz for the ferromagnetic phase, it is easy 
to see that the spin direction of the condensate varies in 
the core of the skyrmion, hence the name.) It is unclear, 
however, whether the spin Skyrmions is stable against 
collapse. 

From the discussion of the SU{?>) phase above, it is 
not a big jump to realize that the ferromagnetic phase 
also supports semilocal strings classified by 7ri(S'^) = Z. 
Furthermore, since the gauge vacuum manifold, which is 
an S^ , is again embedded non-trivially in the full vac- 
uum manifold, we find a fibration of 5*0(3) with fiber S^ 
and base space 0(3)/S'^. Since 7r2 (0(3)75^) = Z, the 
ferromagnetic phase supports global monopoles that sit 
on the ends of the semilocal strings, as does the S'C/(3) 
phase. 

Finally, suppose the SU(3)-broken phase is actually a 
nematic one. The symmetry 0(3) x U{1)em is broken 
down to 0(2) K Z2, and the vortices are classified by 
TTi [0(3) X [/(l)/(0(2) KZ2)], which is not trivial. The vor- 
tex strings in this kind of phase are called "Alice strings" 
and are known to have very unusual properties [42n4B] . 
As one travels around the vortex the direction of the 
condensate changes by tt (not 27r), since the reversed ori- 
entation can be made up by a shift in the phase by tt. 
In the relativistic context, particles moving adiabatically 
around the vortex flip their charge. In our context, they 
flip their spin. 



VI. CONCLUSION 

As we have seen above, the deuteron liquid first dis- 
cussed in Ref. [1] has a rich phase structure. To show 
this, we calculated the one-loop effective potential for 
the deuteron fields, and showed that to leading order. 



n 



C2 >0 


C2 <0 


Ferromagnetic 


Nematic 


(d) =a + zb.ab = 


(d) =ae'",ax b = 


0(3) X U{1)em 


0(3) X U{1)eai 


U{1) 


(9(2) X Za 


SU(3) Phase 


(d) = a + ib 


SU{S) X U{1)em 


U{2) 



FIG. 6: Table of possible phases, their condensates, and cor- 
responding coset spaces as the density n ~ l/l^ increases at 
a fixed temperature T ~ 10"' Ji'. The parameter C2 is the coef- 
ficient of the d^ d^ term in the effective potential, Eq. 
which turns out to be proportional to (^o — A2)/M. 



the theory enjoys an enhanced SU{3) global symmetry. 
We then discussed S'[/(3)-breaking effects, and showed 
that somewhat surprisingly, deuteron-deuteron interac- 
tions mediated by the strong force are the dominant 
source of 5C/(3)-breaking. However, as discussed in 
Sec. |IV| the importance of the 5'C/(3)-breaking effects de- 
pends on the density of the system. At lower densities, 
the system retains the SU{3) symmetry, while at higher 
densities the global symmetry of the deuteron liquid is 
explicitly broken to 0(3). The spontaneous symmetry 
breaking patterns in these two regimes are very different, 
as discussed in Sees. IIIll and FVl 

A brief comment on the relation between the effective 
field theory used here and the one in Ref. [I] is probably 
useful. Our effective theory, valid at distance scales larger 



than the deuteron radius is 'more microscopic' than the 
one in Ref. pj, valid at distances larger than I. As such, 
coefficients that are phenomenological in Ref. [T] can be 
computed in our theory. In addition, we considered a 
more general symmetry breaking pattern (the ferromag- 
netic phase) than addressed in Ref. [1]. 

Our results can serve as foundational work for fur- 
ther exploration of these phases. Let us mention just 
a few possible future directions. In our analysis we as- 
sumed that the condensate is spatially homogeneous, 
but of course it is important to check if a spatially- 
inhomogeneous condensate is also possible. Next, a more 
systematic study of finite temperature effects on the 
deuteron liquid is desirable, as is an explicit computa- 
tion of the condensation temperature. It is also essential 
to better understand the nuclear interactions that break 
the SU{i) symmetry and drive the deuteron liquid into 
either the nematic or ferromagnetic phase. 

Particularly for applications to astrophysics, it would 
be useful to understand the magnetic properties of the 
deuteron liquid. As part of such an investigation, one 
would need to explore under the conditions under which 
the deuteron liquid supports stable vortex strings. More 
generally, it will be interesting to investigate the stability 
and dynamics of the bestiary of topological defects we 
sketched in Sec. |Vl 

Clearly, deuterium at extreme densities turns out to 
be an interesting system deserving of further study. 
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